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Abstract 

Recently, possible indications of line spectral features in the Fermi-LAT pho- 
ton spectrum towards the galactic center have been reported. If the distinct line 
features arise from dark matter (DM) annihilation into 7X (X = 7, Z° or h°), 
the corresponding annihilation cross-section is unnaturally large for typical loop- 
induced radiative processes. On the other hand, it is still too small to be responsible 
for the observed DM relic density. We show that the mechanism of Sommerfeld 
enhancement with scalar force-carrier can provide a simple solution to these puz- 
zles. The possibly large Sommerfeld enhancement of the cross-section for s-wave 
DM annihilation can significantly reduce the required effective couplings between 
DM and charged particles in typical loop diagrams. The DM particles necessarily 
annihilate into scalar force-carriers through tree-level p-wave process, which can 
dominate the total DM annihilation cross-section at freeze out, resulting in the 
correct thermal relic density, but has subdominant contributions to the DM anni- 
hilation today due to velocity suppression. We perform detailed analysis on the 
effects of p-wave Sommerfeld enhancement on freeze out. The results show that 
with the constraints from the thermal relic density, the required effective couplings 
can be reduced by an order of magnitude. 
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1 Introduction 



It has been well-established from observations that dark matter (DM) contributes to 
nearly 23% of the energy budget of the Universe. The leading DM candidates such as 
the weakly interacting massive particles (WIMPs) can interact weakly with the ordinary 
matter and possibly be detected through direct and indirect searches. Monoenergetic 
gamma-ray lines are one of the smoking gun signals of halo DM annihilation, which is 
hard to mimic by astrophysical sources. 

Recently, a possible line spectral feature around 130 GeV or with an additional line 
at 111 GeV in the Fermi-Large- Area- Telescope (Fermi-LAT) photon spectrum in regions 
close to the galactic center (GC) have been reported p~H5]. Indications of similar spectral 
features have been reported in the galactic clusters with modest statistical significance [6] 
and possibly in unassociated gamma-ray point sources [7HTU]. The latest analysis from 
Fermi-LAT collaboration also shows the indication of such line feature at 130 GeV (135 
GeV) with 4.01 (3.35) a local significance using unreprocessed (reprocessed) data [TT] . 
However, no globally significant gamma-ray line features have been established from 
the analysis of Fermi-LAT collaboration. It remains to be confirmed whether these line 
features are indeed from DM annihilation or due to instrumental uncertainties [12] or 
astrophysical backgrounds [T3T[T5] . The upcoming HESS-II experiment can provide an 
independent check on existence of the gamma-ray line features. 

The monoenergetic gamma-ray lines naturally arise from DM annihilation into two- 
body final states jX, where X stands for the Standard Model (SM) neutral parti- 
cles 7, Z° and Higgs boson h° etc.. In this case, the energy of the photon is given 
by Ery = m x [l — m|/(4m^,)], where m x is the DM mass. The observed line sig- 
nals correspond to a thermally averaged cross section multiplied by relative velocity 
(<7 77 v r ci) ~ 1-27 x 10~ 27 cm 3 s _1 for Einasto profile and ~ 2.27 x 10~ 27 cm 3 s _1 for a gen- 
eralized Navarro- Frenk- White (NFW) profile respectively [2] . Although the line spectral 
feature is expected from DM two-body annihilation, the corresponding cross-section and 
non-observation of an excess of accompanying continuum gamma-ray flux may chal- 
lenge the explanation in terms of simple WIMP models [T6HT8] . First, in most DM 
models the DM-photon couplings are generated radiatively through loop level diagrams 
with charged intermediate states EE where F can be SM charged gauge bosons W ± 
and charged fermions / except for the top-quark. For an annihilation cross-section 
(<T 77 v r ci) ~ O(10~ 27 ) cm 3 s _1 , the required DM couplings to mediator particles are in 
general unnaturally large, which may raise the question of perturbativity |16j . Second, 
if the intermediate charged states are kinematically allowed to be the annihilation final 
states such as // and W^W^, the corresponding tree level cross-sections are related to 
that of DM annihilation into 27 as (cr fj^ ww v rc \) / {a~ n v TC \) ~ (7r/a cm ) 2 ps 2 x 10 5 . Such 
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a large cross-section is stringently constrained by the non-observation of any excesses in 
the continuum gamma-ray spectrum and the cosmic-ray antiproton flux [T9H23] . Finally, 
if the XX ~~ FF channels are not opened and XX ~> lX is the main DM annihilation 
channel as suggested by the current observations, the corresponding cross-section is not 
large enough to generate the correct DM thermal relic density which typically requires 
(av re {) p ~ 3 x 10~ 26 cm 3 s _1 for s-wave annihilation. 

Several mechanisms have been proposed to address this problem such as co-annihilation, 
forbidden channel, asymmetric DM [17], resonant annihilation and cascade annihila- 
tion [18] . etc.. Note that many of them introduce degeneracies in the mass of the DM 
particles and the intermediate particles. Specific models in which the signals of XX ~~ >* lX 
and the thermal relic density are not correlated have been considered in Refs. [23H33] . 

In this work, we show that the Sommerfeld enhancement with scalar force-carrier can 
simultaneously explain the large loop level cross-sections and the correct thermal relic 
density without introducing degeneracies in the mass of the DM and mediator particles. 
In the presence of Sommerfeld enhancement, the cross-section for s-wave annihilation 
XX ~~ >" lX today can be greatly enhanced due to the multiple exchange of a light force- 
carrier </>, which reduces the required couplings in the loop diagrams. As the force- carriers 
are light, the DM particles necessarily annihilate into the force-carriers. For being a 
scalar particle, the annihilation XX ~~ >" 00 proceeds through p-wave. The p-wave process 
is also Sommerfeld enhanced, which can dominate the total annihilation cross-section at 
freeze out, lead to the correct thermal relic density, but plays a subdominant role in the 
halo DM annihilation today due to the velocity suppression. The total annihilation in the 
halo can be dominated by 7X channel, which can explain both the reported gamma-ray 
line spectrum and the nonobservation of continuum spectrum. 

This paper is organized as follows. In Sec. EJ we begin with a brief review on the 
main features of the Sommerfeld enhancement, and then focus on the Sommerfeld en- 
hancement of p-wave processes and its impact on the thermal relic density. The effects of 
kinetic decoupling for p-wave processes are also discussed. In Sec. [31 we discuss the max- 
imally allowed Sommerfeld enhancement factor after considering the constraints from 
the thermal relic density for the case with scalar force-carriers. In Sec. HI we apply the 
mechanism to a reference DM model in which the XX ~~ >* %1 proceeds through one-loop 
diagrams and discuss reduction of the required effective couplings between DM particles 
and the charged intermediate states in the loop. The conclusions are given in Sec. [51 
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Figure 1: (Left) Feynman diagram of DM annihilation into 27 with multiple-exchange 
of force-carriers which leads to the Sommerfeld enhancement; (Middle) diagram of t- 
channel DM annihilation into force-carrier particles; (Right) diagram of s-channel DM 
particle annihilation into force-carrier particles through <p 3 type interactions. 

2 Sommerfeld enhancements of £>-wave DM annihi- 
lation and thermal relic density 

The Sommerfeld enhancement of annihilation cross-section occurs when the annihilating 
particles self- interact through a long-range attractive potential V(r) at low velocities 
[34] . In this scenario, the short-distance DM annihilation cross-section can be greatly 
enhanced due to the distortion of the wave functions of annihilating particles at origin 
[351438] . The attractive potential may originate from multiple-exchange of light force 
carrier particles between the annihilating DM particles as shown in Fig. [U The nature 
of Sommerfeld enhancement have been extensively studied (see, e.g., Refs. [39144*8] ) in 
light of the cosmic-ray positron/electron excesses observed by PAMELA [19], ATIC [50J, 
and Fermi-LAT [5T] etc.. 

In the previous studies, the force-carrier is often assumed to be a vector boson, as 
it can be naturally light [39]. The maximally allowed enhancement factor turned out 
to be stringently constrained by the thermal relic density due to the additional s-wave 
annihilation XX ~~ >* 00 i n t ne presence of the light force-carrier , which may challenge 
the Sommerfeld enhancement as an explanation for those excesses [4Tll4*8l l52] . Note that 
the cosmic-ray lepton excesses may have astrophysical origins and may not be related to 
the halo DM annihilation. In this work, we shall consider the Sommerfeld enhancement 
as a mechanism to simultaneously account for both the possible gamma-ray line signals 
and the correct thermal relic density. For this purpose, we shall focus on the case of 
scalar force- carriers. Although in both the vector and scalar cases the induced long-range 
attractive potential is of the same Yukawa type. For a scalar force-carrier, the XX ~^ 00 is 
dominated by p-wave processes, which significantly modifies the relic density constraints 
due to the different velocity dependencies. The two-body wave function \l/(r) of the 
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annihilating DM particles satisfies the non-relativistic Schrodinger equation 

1 



m 



-V 2 ^(r) + V (r) (r) = m x f * ( r ) 



where v = v re \/2 is the velocity of DM particle in the center-of-mass frame and v Te \ is 
the relative velocity of the annihilating DM particles. After an expansion over angular 
momentum £, namely, \I/(r, 9) = '^2 £ Pe(cos8)xe{r)/r, the Schrodinger equation for radial 
wave function xe can be written as: 



d\e (t) 
dt 2 



£(£+1) + V(t) 



t 2 



m x v A 



Xe {t) + xe (t) = 0, 



(2) 



where t = m x vr. The above Schrodinger equation can be solved with the following 
boundary conditions pTljl^2"] 



lim xe (t) = t l+1 and lim xt (t) — > C sin I t — 



t-to 



t—too 



£% 



8, 



(3) 



where 5g is the phase shift and C is a normalization constant. With the above boundary 
conditions, the Sommerfeld enhancement factor Sg of the annihilation cross-section is 
given by [39] 



S, 



lim 

t->o 



Xe (t) 
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~(2£+ l)!f 


xf (t) 




C 



(4) 



where xf\t) * s the wave function in the free- motion case without a potential. The 
exchange of massive vector or scalar particles between the DM particles results in an 
attractive Yukawa potential 

V(r) = , (5) 



where a = g 2 /{An) is the coupling strength of XX'P or X7 M X</V type of interactions. The 
nature of the Sommerfeld enhancement depends on the two variables 



s v = — and es 

a 



am 



(6) 



x 



In the limit of < ej, the Yukawa potential can be well approximated by a Coulomb- 
type potential. The corresponding Schrodinger equation can be solved analytically for 
arbitrary angular momentum and the enhancement factors are 



S 



Col 



r-Col 
°0 



1 



exp (-7r/e„)' 
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71=1 



An 2 e 2 



(for £ = 0), 



(for £^0). 



(7) 
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For small e v /ix <C 1, the enhancement factors can be approximated by Sf o1 ~ 2ir / ((2e v ) 2£+1 (£\) 2 ). 
Therefore, at low velocities the s- and p-wave Sommerfeld enhancement factor scales as 
1/v and (1/f 3 ) respectively. 

In the case where is non-negligible, the 1/v behavior of s-wave cross-section breaks 
down. Through approximating the Yukawa potential by the Hulthen potential, the s- 
wave Sommerfeld enhancement factor can be estimated as [^2|H3] 



For ^> e v , namely, the deBroglie wavelength of incoming particles is longer than the 
range of the interaction, the s-wave Sommerfeld enhancement saturate with Sq ~ 12/e^. 
But for some particular values of ~ 6/(n 2 n 2 ), (n = 1,2,3,...) for which the DM 
can form zero-energy bound states, there exists additional resonant enhancements which 
scale as 1/v 2 . The resonant enhancement is eventually cut off by the finite width of the 
resonance [39] . 

For p-wave Sommerfeld enhancement factor Si, there is no analytic expression avail- 
able. The p-wave case has been investigated in Refs. [^0l - H2] without considering its 
effects on the freeze out and thermal relic density. In this work, we shall focus on these 
effects as they are important in determining the maximally allowed Sommerfeld enhance- 
ment. We first numerically solve the Eq. (j2J) with the boundary conditions in Eq. (J3J) 
and illustrate the dependence of Si on the two variables e v and in Fig. |2j For small 
e 4> ~ 10~ 3 , the value of Si scales as 1/v 3 as expected from the Coulomb limit. For larger 

in the range 10~ 3 — 10 _1 , it has resonant behavior which is similar to the s-wave case. 
For even larger > 10 _1 the enhancement saturates. Note that the generic p-wave an- 
nihilation cross-section before including the Sommerfeld enhancement is proportional to 
v 2 , thus the velocity dependence of the total Sommerfeld enhanced p-wave annihilation 
cross-section should be proportional to Sie 2 . In the right panel of Fig. [21 the contours 
of Si el in the (e v ,e < j > ) plane are shown. In the region where < 10~ 3 , it scales as 1/v 
instead of l/t> 3 . In the resonance region 1CT 3 < < 1CT 1 , the velocity dependence 
of Sie„ is not significant. In the saturation region > 10 _1 , S^ 2 , scales as v 2 , the 
total cross-section decreases rapidly towards low velocities. Thus the main difference 
from the s-wave case is that the total p-wave annihilation cross-section can be either 
velocity-suppressed or velocity-enhanced, depending on the size of e^. 

An important quantity directly related to the thermal relic density is thermally av- 
eraged annihilation cross-section. The generic DM annihilation cross-section times the 
relative velocity before including the Sommerfeld enhancement has the form (av rc i)o = 
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Figure 2: (Left) Contours of p-wave Sommerfeld enhancement S\ with Yukawa potential 
in (e„,e^) plane. (Right) the same but for Sie^ which is relevant to the total p-wave 
annihilation cross-section. 

a + k^ cl + 0(vf cl ), where a and b are coefficients correspond to the s- and p-wave contri- 
butions which are velocity-independent. After including the Sommerfeld enhancement, 
the thermally averaged cross-section at a temperature T can be written as 



where the thermal average of a quantity X(v re i) in the non-relativistic limit is given by 



with x = m x /T. The thermally averaged annihilation cross-section is a function of x 
and depends on the parameters a and m^. 

The time evolution of the DM number density n x is governed by the Boltzmann 
equation 



where y( eq ) = n x / 's is the (equilibrium) number density rescaled by entropy density s, 
mpi ~ 1.22 x 10 19 GeV is the Planck mass scale, g^s and g* are the effective relativistic 
degrees of freedom for entropy and energy density respectively. The decoupling temper- 
ature Xf is defined as the temperature at which the dark matter particles start to depart 




(9) 




(10) 
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from the thermal equilibrium, and the density Y is related to the equilibrium density 
yeq \jy y (xf) = (1 + c)Y cq (xf), where c is a constant of order unity. The value of Xf is 
approximately given by [53] 




with g x the degrees of freedom of dark matter particle. In the absence of Sommerfeld 
enhancement (<7f rc i) = a + 6b/xf, taking c ~ 1(2) for s(p)-wave annihilation leads to 
good fits to the numerical solutions of the Boltzmann equation. The DM number density 
in the present-day can be obtained by integrating Eq. ( 112]) with respect to x in the region 
Xf < x < x s , where x s corresponds to the temperature at which the DM annihilation rate 
is insignificant compared with that of the expansion of the Universe and Y{x) becomes 
stable. The value of Y(x s ) can be written as 



In performing the integration over x, as {<Jv ve \) depends on temperature, one needs to 
take into account the effects of kinetic decoupling. When the DM particles are in both 
chemical and kinetic equilibrium with the radiation background, the temperature of the 
DM particles tracks that of the background, i.e., T x = T or x x = m x /T x = x. After 
dropping out of chemical equilibrium, the DM particles can still remain in kinetic equi- 
librium with the radiation background through scattering off SM relativistic particles 
which are in thermal equilibrium with the radiation background. At some temperature 
Tkd, when the rate of the scattering cannot compete with that of the expansion of the 
Universe, the DM particles start to decouple from kinetic equilibrium. After the kinetic 
decoupling, T x drops quickly with the scale factor a as a~ 2 instead of a -1 , the temper- 
atures of the DM particles and the radiation background are approximately related by 
T x = T 2 /T kd , or x x = x 2 /xkd [M]- Thus the integration from Xf to x s needs to be 
separated into two parts, from Xf to x^d and from x^d to x s . For the second part of the 
integration, one should use x x instead of x. Previous analysis have shown that including 
the effect of kinetic decoupling leads to a significant reduction of relic density in the case 
of s-wave annihilation with Sommerfeld enhancement [4"5ll4"6| l4"8]. The value of T^d is 
model-dependent, for instance, in supersymmetric models ps 10~ 3 — 10 _1 Tj [55]. In 
this work, we shall take the value of free parameter. Finally, after freeze out, 

the relic abundance of DM particles is given by 




x f ps ln[0.038c(c + 2)m Pl m x (g x g^ 1/2 )(crv Tcl )} 

-- lnln[0.038c(c + 2)m P im x (g x g~ 1/2 )(crv TC i)}, 



(13) 




(14) 




(15) 



which is to be compared with the observed value Qh 2 = 0.113 ± 0.004 |56j . 

We solve the Boltzmann equation Eq. (fl~2l) numerically including the effects of Som- 
merfeld enhancement and kinetic decoupling for both s-wave and p-wave annihilation. 
In order to facilitate the comparison, we take a = 2.2 x 10 -26 cm 3 s _1 and b = 1.7 x 
10 -25 cm 3 s _1 such that the final DM relic abundance Qh 2 ~ 0.113 is the same for both 
s-wave and p-wave in the absence of Sommerfeld enhancement and kinetic decoupling. 
We fix m x = 130 GeV and = 0.25 GeV and consider two different values of a = 0.1 
and 0.01, and two kinetic decoupling temperatures Tf = 2Tkd and lOTkd- The evolutions 
of Y(x)-Y e i{ X^j clS cl function of x are shown in Fig. [3j For both s- and p-wave annihila- 
tion, the inclusion of Sommerfeld enhancement results in reduction of final thermal relic 
density by a factor of (9(1). For a = 0.1 the relic density is reduced by a factor of ~ 3 for 
s-wave annihilation while ~ 2 for p-wave case. At high temperatures x ~ 20, the p-wave 
annihilation cross-section still decreases with temperature in the presence of Sommerfeld 
enhancement, which leads to earlier decoupling from the thermal background and larger 
relic density. As it can be seen in Fig. [31 the effects of kinetic decoupling are significantly 
different between s- and p-wave case. For Tf = 2Tkd and lOTkd, the kinetic decoupling 
leads to further reduction of relic density by ~ 50%(20%) in the s-wave case. However, 
for p-wave case, the reduction is almost invisible for a = 0.1 as shown in Fig. [3j For 
smaller a = 0.01, it even leads to a slight enhancement of the relic density. This is 
because the kinetic decoupling makes the DM particle freeze out more quickly in the 
p-wave case as the annihilation cross-section decreases with temperature. 



3 Sommerfeld enhancement of XX ~^ 1^ cross-sections 

With the main features of p-wave Sommerfeld enhancement given in the previous section, 
we are ready to discuss the contribution to the DM thermal relic density from the process 
of DM particles annihilating into scalar force-carriers through p-wave. Since the force- 
carrier is much lighter than the DM particle, i.e., <^ m x , the DM particles necessarily 
annihilate into the force-carriers, which contributes to a DM annihilation channel in 
addition to jX, and can be the dominant contribution to the total DM annihilation 
cross-section. The process itself is also Sommerfeld-enhanced, which complicates the 
calculations. 

Before including the effects of Sommerfeld enhancements, the total DM annihilation 
cross-section (atot^rci)o can be written as 

Kot^rcOo = 0<^rcl)o + (^ 7 X^rel)o- (16) 

The DM particles can always pair-annihilate into <fxft through t-channel x-exchange. In 
addition, if there exists non- negligible cubic and quartic self-interactions between the 
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s-wave, m =130GeV, m =0.25GeV,a=0.1 

z ♦ 



p-wave, m =130GeV, m =0.25GeV, a=0.1 

x ♦ 




Figure 3: Effects of Sommerfeld enhancement and kinetic decoupling on the temperature 
evolution of the DM number density for the case of s-wave (left) and p-wave (right) with 
a =0.1 and 0.01, see text for explanations. 



force- carriers of the form — /i0 3 /3! — A0 4 /4!, two-body or three-body s-channel annihila- 
tion may occur. The three-body annihilation is highly suppressed by small phase-space 
and is neglected in this work. In the case of scalar force- carrier, both the t- and s-channel 
two-body annihilation into cfxp corresponding to the two diagrams in Fig. [1] are p-wave 
processes. In the limit of <C m x , the total annihilation cross-section of XX ~~ >* 00 
from the calculation of the two diagrams in Fig. [I] is independent of m^, and is given by 

, . 37ra 2 / 5 . 1 2 

fo^)°= 8^ 1_ i8* + 48* (17) 

with £ = /i/ (2m x y/ an) . The first term in the right-hand-side of Eq. ( 1TT1) comes from 
the t-channel contribution. The term proportional to £ 2 corresponds to the squared 
amplitude of s-channel annihilation, and the term proportional to £ corresponds to the 
interference between s- and t-channel diagrams. We assume that the other annihilation 
channel XX ~~ lX is an s-wave process such that there is no explicit velocity dependence 
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in (avyjf^reiV After including the Sommerfeld enhancement, the total thermally averaged 
annihilation cross-section is 



(Otot^rel) 



(c 7 X^rel)o(<So) + ((o"^ re i)o<S'i). 



(18) 



The Sommerfeld enhancement factors S^i depend on parameters a and m^, and the 
cross-section (o"^w re i)o is also a function of a. Therefore, the requirement of reproducing 
the correct thermal relic density determines the size of a for a given (cr^x v rc i)o and 
force-carrier mass m^. Using the value of a constrained by the thermal relic density, 
one can calculate the allowed Sommerfeld enhancement factor (So) for (cr 7 x^rci)o at l° w 
temperatures, which is needed to predict the cross-section of XX ~^ lX in the present 
day 
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Figure 4: The value of a allowed by the DM thermal relic density as a function of 
((X-yX^reiV (Left) the cases without kinetic decoupling. Three curves correspond to 
H — (solid), m x /2 (dashed), and m x (dot-dashed), respectively (Right) the case 
with kinetic decoupling. Three curves correspond to without kinetic decoupling (solid), 
Tkd = Tf/10 (dashed) and Tf/2 (dot-dashed), respectively with /j, fixed at m x . 



In the left panel of Fig. HI we show the allowed values of a as a function of (cx^ u re i)o 
for a fixed = 0.25 GeV with three different choices of cubic coupling [i. For sufficiently 
small (a-yx^reOo <^ 1CT 26 cm 3 s _1 , namely, the total annihilation cross-section is domi- 
nated by (a^v ie \)o at freeze out, the value of a is found to be insensitive to (cr 7 xv r ei)o- 
From the figure, one obtains a ~ 0.014, 0.016, and 0.019 for \x = 0, m x /2, and m x , re- 
spectively. When ((T 7 x^rei)o grows and approaches ~ 10~ 26 cm 3 s -1 , the value of (cr^v^o 
needs to be suppressed in order to give the correct relic density, which results in a sig- 
nificant reduction of a. As it can be seen from the figure, at (<7 7 x^rci)o ~ 10 -26 cm 3 s~ 1 , 
the allowed values of a is reduced by ~ 50%. 
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The effects of kinetic decoupling are shown in the right panel of Fig. H]for two different 
decoupling temperatures Ty = 2T kd and 10T kd . For small (cr 7 xf r ei)o ^ 10~ 26 cm 3 s _1 , 
the inclusion of kinetic decoupling only leads to a slight increase of ct, which is expected 
as in this region the p-wave annihilation dominates, the kinetic decoupling makes DM 
particle freeze out more quickly, which leads to a higher relic density as discussed in the 
previous section. For large (o" 7 xf rc i)o ~ 10~ 26 cm 3 s _1 , the effect of kinetic decoupling 
is sizeable, as in this region, the s-wave annihilation becomes important, which has a 
stronger dependence on the kinetic decoupling than that of p-wave case. This leads to a 
significant decrease of a. Especially in the vicinity of resonances, due to the additional 
s-wave resonant Sommerfeld enhancement, there could be a sudden reduction on the size 
of a, which can be seen at (f7 7 x u rei)o ~ 6 x 10~ 27 cm 3 s~ 1 in Fig. HI 

After constraining the allowed values of a, the s-wave Sommerfeld enhancements fac- 
tor (So) at low temperatures can be calculated straightforwardly. The velocity-averaged 
Sommerfeld enhancement of the halo DM annihilation cross-section in the present-day 
is defined as 

1 [2 f VcBC -^ cl 
(So)now = T7-W- / Sb e 2 *o v 2 reX dv rc i , (19) 

where Vq is the DM velocity dispersion and v esc is the DM escape velocity, respectively. 
N = erf (z/y/2) ~ (2/vr) 1/2 ze~ z2 / 2 is a normalization constant with z = v csc /vq. Both v esc 
and Vq depend on the distance r from the GC. In the vicinity of the Sun, r = r « 8.5 
kpc, f eS c( r 0) ~ 525 km s -1 , and v (r®) ~ 210 km s -1 . In the left panel of Fig. 
we show the dependences of (So)now on the coupling a and the velocity dispersion t> . 
For the allowed values a = 0.0136, 0.0159, and 0.0186 corresponding to fi = 0, m x /2, 
and m x , the enhancement factors are (So) now ~ 170, 70, and 50, respectively. Thus 
the Sommerfeld enhancement factor can reach (9(100) with the constraints from DM 
thermal relic density, which is larger than the case where the force-carrier is a vector 
boson PUSH]. 

We are now in the position to discuss the line spectral shape in the photon spec- 
trum recently observed by the Fermi-LAT collaboration. The observation, if inter- 
preted as DM annihilation into 77, corresponds to an s-wave velocity-averaged cross- 
section (<7 77 f rc i) 

now °° 10 cm s . All of our previous results on the annihilation 
XX ~^ lX can be directly applied to the case of XX ~^ 77; as it is a special case 
of jX. In the right panel of Fig. [51 the relation between (cr 77 w re i) now and (cr 77 u rel ) 
is shown. One sees that in order to reproduce the observed signal corresponding to 
(o" 77 f rc i) now ~ 1.27 x 10~ 27 cm 3 s _1 for the Einasto profile [2], the required cross-sections 
before the Sommerfeld enhancement can be quite small (cr 77 w re i)o=7.2 x 10~ 30 cm 3 s _1 , 
1.8 x 10~ 29 cm 3 s _1 , and 2.5 x 10~ 29 em's" 1 , for \i — 0, m x /2, and m x , respectively. 
The cross-section of XX ~^ 77 ma Y De further enhanced due to the possible lower DM 
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Figure 5: (Left) (So) 

now 

function of coupling strength a for the velocity dispersion 
vq = 210 km s _1 , 100 km s _1 , and 50 km s _1 respectively. (Right) the relation between 
(o" 77 f rc i)now and (a 77 f rc i)o for three different values of p. The gray band denotes the fit 
cross-section (cr 77 u) = (1.27 ± 0.32) x 10~ 27 cm 3 s _1 in the case of Einasto profile [2]. 
Other parameters are fixed at m x = 130 GeV and = 0.25 GeV. 



velocity dispersion near the GC than that in the solar neighborhood. The dependence 
of Sommerfeld enhancement on v (r) at GC has been discussed in Refs. [571158]. The 
N-body simulations suggest that Vo(r) is related to the DM density profile p(r) through 
a relation vo(r) 3 / p(r) m r x [59] . From pure DM simulations, the power index is found to 
be x - 1-9 - 2.0 [60]. For the NFW profile p(r) oc (r/r s )~ a (l + r/r s y 3+a with a = 1.0 
and r s = 20 kpc, Vo(r) scales as r x ~ a . Therefore it decreases with the distance r. This is 
also true for the Einasto profile p(r) oc exp[— (2/a)(r/r s ) a ], if one takes the typical value 
of a = 0.17. Thus a larger s-wave Sommerfeld enhancement at GC can be expected. 
However, if the Sommerfeld enhancement saturates at some velocity smaller than vo(r & ), 
then it is insensitive to the choice of v (r). The additional enhancement at GC in this 
case can be simply estimated as S (v = 0)/S (v (r )) [61]. Note that the simulation 
results can be modified significantly after including baryons [62] . 

Although at freeze out the p-wave cross-section (a^v re i) can be a few order of magni- 
tudes larger than the s-wave cross-section (cr 77 t> rc i), at lower temperatures it is possible 
that (cr^frci) becomes comparable with or even smaller than (cr 77 t> rc i), which is due to 
the dramatic difference in velocity-dependencies between s- and p-wave processes in the 
presence of Sommerfeld enhancement. This possibility is shown in Fig.[6j At the temper- 
ature of thermal decoupling x ~ 25, one sees that (cr^v re ]) « 1 x 10 -25 cm 3 s _1 which is 
about four order of magnitudes larger than (cr 77 v rc i) • When the temperature goes down, 
in general the value of (a^v^i) decreases due to the velocity-suppression. The p-wave 
Sommerfeld enhancement of (a^Urei) can be significant in the range 10 4 < x < 10 6 when 
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it is near a resonance region. But for other values of x which are off the resonance, 
(0>0frei) decrease rapidly. On the other hand, the size of (cr 77 f re i) increases monotoni- 
cally towards larger x due to the s-wave Sommerfeld enhancement. At x ~ 25, the value 
of ((TjjVrei) is ~ 10~ 29 cm 3 s _1 , but at x ~ 4 x 10 6 which corresponds to vo = 210 km s _1 , 
it reaches ~ 10~ 27 cm 3 s _1 and becomes comparable with (o^iVei) for ji = m x /2 and 
m x . For /i = case, it can even dominate over {(J^v TC \). Therefore, it is possible 
that XX ~~ >* 77 can be the main DM annihilation channel today, and is responsible for 
the observed gamma-ray line at the GC. In this work, we consider a light scalar with 
77i0 ~ 0.25 GeV, such that <fi can decay into fi + fi~ and e + e~. But the final states of 
7r°7r° and t + t - are kinematically forbidden, which will suppress the generation of large 
continuum gamma-ray flux from the annihilation XX ~^ 00 followed by the decay of 
into these final states. 

Note that such a large s-wave enhancement of (9(100) is still consistent with the 
limits derived from BBN and CMB data (63], as the total DM annihilation cross-section 
after the Sommerfeld enhancement is still an order of magnitude smaller than the typical 
WIMP annihilation cross-section. 




x=m fT 

x 



Figure 6: Temperature evolution of (o"^f rc i) and (cr 77 f re i) for three different values of 
\i = 0, m x /2, and m x , respectively. 



4 Sommerfeld enhancement in a simple model 

We have shown in the previous section that the Sommerfeld-enhanced cross-section 
((TyyVrei) ~ 10~ 27 cm 3 s _1 today can be consistent with the DM thermal relic density, 
due to the p-wave annihilation of DM particles into the force-carriers. In this section, 
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Figure 7: Feynman diagram for XX ~> 77 through loop diagrams with an intermediate 
pseudoscalar A and a charged fermion / in the loop. 



we discuss another advantage of the Sommerfeld enhancement, namely, it can also solve 
the problem of unnaturally large couplings required by many DM models motivated to 
explain the putative gamma-ray line. As the DM particles carry no electric charge, it 
is often assumed that the annihilation of DM particles into 77 proceeds only through 
one-loop diagrams which involve massive charged particles running in the loop ( for ex- 
ceptions, see e.g. Refs. [HUE!] ). Due to the loop suppression, in general the effective 
DM couplings to the charged particles in the loop have to be above order unity, in or- 
der to give the cross-section suggested by Fermi-LAT (see e.g. [IE]), which may raise 
the issue of perturbativity. Invoking the mechanism of Sommerfeld enhancements can 
significantly reduce the required couplings to the perturbative region. 

For a concrete illustration, we consider a realization of Sommerfeld enhancement in 
a simple reference model. In this model, the DM particle x is assumed to be a Majorana 
fermion. Other particles in the models are: a light scalar force-carrier 0, a heavy pseu- 
doscalar mediator A , and a heavy charged Dirac fermion / with electromagnetic charge 
number Qf and color number Cf. The relevant interactions in the model are given by 
the following Lagrangian 



-f X0X - *yX75X^° - ig f h 5 fA°. 



(20) 



In order to avoid constraints from the nonobservation of continuum photon spectrum, 
we consider the case where m^m^ 3> m x , such that XX ~~ >* ff,A°A° are kinematically 
forbidden, and XX can only annihilate into 77 through loop process. The corresponding 
Feynman diagram is shown in Fig. [7J The annihilation cross-section for this process is 
given by [T7] 



(cr 



77 



Vre\) 



a: 



m 



47T 3 (s - m 2 A ) 2 + u? A Y\ 



-1 4 



arctan 



^Jm 2 /m 



(21) 



where s is the center-of-mass total energy. The total width of the pseudoscalar A 
receives contribution from the decay channels A — > XX-, J I '■> an d 77 with the partial 
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widths 



ghn A \ 4m 2 g 2 f m A / 4m 2 , 

xx 16tt V m 2 , ' " 8tt V m 2 , ' V ; 



and 



■in 



A a2 em9}QjC 2 f 



r - = 2 5 7:q (23) 



where the function Ay 2 is defined as A^, 2 (r) = 2r 1 /(r) with 



^ 2 io uciuicu cw vi-^ 

{arcsin 2 -y/r (for r < 1) 

-i In — v , - Z7T (for T > 1) 

We assume that in the dark sector the parity symmetry is conserved such that the decay 
A — >• (j)(p is forbidden. Since A is a pseudoscalar, the annihilation XX ~~ >* 77 proceeds 
through s-wave. For heavy rriA, mf 3> the srze °f the annihilation cross-section can 
be estimated as 

(rr „ \ in-^3,-1 (9x9f\ 2 ( m x V f 500 GeV \ 4 / 500 GeV \ 2 

One sees that a very large product of the couplings y/g x 9f ~ 10 is required, which is 
unnatural and can invalidate the perturbative calculations. This is a common problem 
in many dark matter models in which the DM annihilation into 77 through one-loop 
diagrams. 

One way to enhance the cross-section is to assume that the annihilation is near a 
resonance which occurs if m A « 2m x + 5 with 5 <C m A . In this case the cross-section 
can be enhanced by a factor of m 2 A / (4<5 2 + T A ). If small couplings g x , gf ~ 1 is required, 
the relative mass difference should be around 8/m A ~ 2% for « r xx . Another 
modest enhancement may arise from the case where m x ps mf such that the function 
arctan[(m 2 /m 2 — l) -1 ^ 2 ] reaches its maximal value n/2. Compared with the case where 
mf « 5m x the enhancement of the value of the arctan is ~ 8. In addition, in order to 
have the correct relic density, it is required that 8 must be positive [IS] . 

With the presence of Sommerfeld enhancements, the required couplings can be re- 
duced significantly without introducing any mass degeneracies. In Fig. [HI we show the 
required couplings which can reproduce the observed (cr 77 f rc i) = 1.27 x 10~ 27 cm 3 s~ 1 . 
Compared with the case without Sommerfeld enhancements, the required products g x gj 
can be reduced by an order of magnitude. For a wide range of the pseudoscalar mass 
100 < m A < 450 GeV, the required couplings are smaller than unity. If the Sommerfeld 
enhancements are absent, only a narrow region close to 2m x can be consistent with a 
small g x g f . 
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u=m 12 
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fi=m 
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Figure 8: The required value of g x gj to account for the gamma-ray line signals with 
cross-section (cr 77 u re i) = 1.27 x lCr 27 cm 3 s -1 [2J. Three cases with Sommerfeld enhance- 
ment correspond to ji = (solid), m x /2 (dashed), and m x (dot-dashed), respectively. 
The dotted line denotes the case without Sommerfeld enhancement. Other parameters 
are fixed at Qf — Cf — 1, g x = gj, and rrif = 300 GeV. 



In this simple model, the force-carrier does not appear in the leading one-loop 
diagram for XX ~^ Hi which make it straightforward to factorize the short- and long- 
distance contributions. In general, the force-carrier can play important roles in both 
the one-loop diagram and the long-distance Sommerfeld enhancements. For instance, 
in the case where the DM particle is a Majorana fermion belonging to the adjoint pre- 
sentation of SU{2) W gauge group (e.g Wino neutralino), the SU{2) W gauge bosons W ± 
appear in both the box-diagrams for DM annihilation into 77 and the long-distance 
ladder diagrams as a force-carrier which results in off-diagonal Sommerfeld enhance- 
ments [35H38| B4,66j. A fully consistent analysis of incorporating these two contribution 
simultaneously can be found in Ref. [67] . 

5 Conclusions 

In summary, the recently reported possible indications of line spectral features in the 
Fermi-LAT photon spectrum in regions close to the galactic center can be related to the 
signals of halo DM annihilation. However, the corresponding annihilation cross-section 
of (cr 77 t' r ci) ~ O(10~ 27 ) cm 3 s _1 is too large for typical loop-induced radiative processes, 
while on the other hand too small to be responsible for the observed DM relic density 
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which is typically 3 x 10~ 26 cm 3 s _1 . We have pointed out that the Sommerfeld en- 
hancement with scalar force-carrier can simultaneously explain those features. In this 
mechanism, the cross-section for DM annihilation into jX in the galactic center today 
can be greatly enhanced due to the attractive forces between the DM particles, which is 
induced by the multiple-exchange of the force- carriers. The additional p-wave annihila- 
tion into the force-carriers can dominate the total DM annihilation cross-section at freeze 
out and set the correct thermal relic density, but can have subdominant contributions 
to the DM annihilation today due to the velocity suppression. We have performed an 
detailed analysis on the allowed Sommerfeld enhancement under the constraints from 
the thermal relic density determined by p-wave annihilation processes. The results show 
that the Sommerfeld enhancement factor can reach (9(100). In a simple reference model 
with XX ~~ 27 occurring at one loop, we have shown that the required DM effective 
couplings to the intermediate particles in the loop can be reduced by an order of mag- 
nitude and below unity, which keeps the perturbativity of the model. Compared with 
some other mechanisms for increasing the DM annihilation cross-section at loop level, 
the Sommerfeld enhancement does not require any degeneracies in the masses of DM 
particles and the intermediate particles in the loop diagrams. 
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